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High rates of recession in arc jet shear tests of Phenolic Impregnated Carbon Ablator 
(PICA) inspired a series of tests and analysis on FiberForm (a carbon preform used in the 
fabrication of PICA). Arc jet tests were performed on FiberForm in both air and pure 
nitrogen for stagnation and shear configurations. The nitrogen tests showed little or no 
recession, while the air tests of FiberForm showed recession rates similar to that of PICA 
(when adjusted for the difference in density). While mechanical erosion can not be ruled 
out, this is the first step in doing so. Analysis using a carbon oxidation boundary condition 
within DPLR was used to predict the recession rate of FiberForm. The analysis indicates 
that much of the anomalous recession behavior seen in shear tests may simply be an artifact 
of the non-flight like test configuration (copper upstream of the test article) a result of dis- 
similar enthalpy and oxygen concentration profiles on the copper. Shape change effects 
were also investigated and shown to be relatively small. 


Nomenclature 


C H = Stanton number for heat transfer, qj(p e u e (h e h w )) 

C M = Stanton number for mass transfer, m" /(( p e U e )(c 5 — C s )) 

c s = mass fraction of species s 

h = enthalpy (MJ/kg) 

M = Mach number 

p = pressure (kPa) 

q = heating rate 


I. Introduction 

C arbon ablators are proposed for numerous planetary entry missions, such as Mars Science Laboratory 
(MSL), Crew Exploration Vehicle (CEV), and Venus missions. Carbon ablators typically consist of a 
carbon fiber matrix that is infiltrated with a phenolic polymer such as is the case with Phenolic Impregnated Carbon 
Ablator (PICA) 1 or Carbon Phenolic 2 . 

In order to qualify candidate TPS materials for flight, testing is performed over a range of conditions in an arc jet 
in both stagnation and shear (often in a wedge test fixture). In the case of PICA the recession rate is often greater in 
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Figure 1 Measured Recession of PICA in a shear 
environment relative to the FIAT prediction. 


shear than it is in a comparable stagnation test 3 . The 
reason for this difference has been unknown. 

Consequently, heat shield designers 4 have been forced 
to add margin to account for the higher than predicted 
rate of recession seen in shear. Figure 1 shows the 
measured recession rate normalized by a recession 
rate prediction by the Fully Implicit Ablation and 
Thermal Response Program (FIAT) 5 , for both 
stagnation and shear tests. The ratio of measured to 
FIAT predicted recession for stagnation is 
approximately one, while for shear tests it is about 1.5 
on average, and in some cases the ratio is 2 (and as 
much as 2.5 when uncertainty is included). 

Several theories have been proposed for the 
difference between shear tests and stagnation tests; 

1) Mechanical erosion 

2) Shape change effects 

3) Oxygen transport to the wall is not proportional to heat transfer (as assumed in FIAT). 

It is worth further expanding on the 3 rd theory that oxygen transport is not proportional to heat transfer. Shear tests 
often rely on a water-cooled copper fixture to hold the test specimen and protect it from the extreme heating 
associated with stagnation at the nose (typically 10 times that of the test article). In the case of a wedge, the flow 
stagnates on the copper nose where the heating can 
reach 3000 W/cm 2 due to the small nose radius. The 
arc-heated air travels for some distance over the 
copper fixture developing a thermal boundary layer 
that is in equilibrium with the cold copper wall and a 
thermal gradient normal to the wall. In this same 
upstream copper portion of the test fixture, oxygen is 
in contact with the wall, but is not being “consumed” 

(oxygen reacting with carbon) and therefore there is 
no oxygen gradient normal to the wall (where oxygen 
can be either O or 0 2 ). 

Consequently, the flow over the test article is 
oxygen rich (relative to configurations that have 
carbon leading edges). In contrast, the enthalpy is 
depleted near the wall (relative to the free stream), due 
to contact with the water-cooled (cold) copper fixture. 

This lack of similarity between the enthalpy and 
oxygen gradients, generated on the copper fixture 
upstream, persists down stream (see Fig 2), and leads 
to a mismatch between heat transfer coefficient and 
oxygen mass transport coefficient. FIAT assumes the 
two coefficients are equal. This effect could be 
responsible for the high rate of recession seen in the 
experiments that is not accurately predicted by FIAT. 

Coupling between the flow field and surface reactions 
is useful for studying non-canonical boundary layer 
flows, such as occurred in the shear tests. The reason 
for this will become apparent in the following section. 

This paper describes an effort to understand the 
relatively high rate of recession seen on carbon 
ablators in a shear environment, by exploring each of 
the aforementioned theories. 
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II. Arc Jet Tests of FiberForm 

In the interest of testing theories as to why the material response codes are under-predicting the recession in 
shear tests, several arc jet tests were performed on FiberForm©. FiberForm is a lightweight carbon structure which 
is the carbon preform into which phenolic is impregnated in the course of making Phenolic Impregnated Carbon 
Ablator (PICA). Since FiberForm is the carbon backbone of PICA it will likely have similar properties to that of the 
char layer of PICA (properties such as conductivity, density, tensile strength). Stackpoole et al. 6 measured the PICA 
char layer and found that the density of the outermost layer of char was equal to that of FiberForm. 

The main purpose of the phenolic is to 1) provide a sink for energy to slow the flow of heat (via conduction) to 
the bond-line and substructure of the vehicle and 2) plug the pores of the FiberForm and inhibit the flow of hot 
gasses through the FiberForm. Because FiberForm does not contain phenolic, it is much simpler to model than 
PICA - it has fewer chemical reactions and no pyrolysis gasses. Never the less the recession behavior of FiberForm 
in arc jet tests is similar to that of PICA because most of the recession is due to carbon oxidation. FiberForm 
recesses at a 50% greater rate than PICA mostly due to its 33% lower density. 

Many arc jet tests were performed on FiberForm and several are briefly described in this section. The data is 
proving useful for understanding the high recession rates seen on PICA in shear tests (described later in this paper). 


A. Stagnation Tests 

One such test was 
performed in the 20MW 
Aerodynamic Heating 

Facility (AHF) 6 arc jet at 
NASA Ames Research 
Center using nitrogen as the 
test gas in one case and air as 
a test gas in another (see Fig. 
3). These tests were 
performed on 4” IsoQ shaped 
stagnation models in the 7” 
nozzle of the AHF at near 
heater 


maximum 


arc 


FiberForm 
AHF-281-F-03 
120 sec 


FiberForm 
AHF-281-F-02 
60 sec 




a) FiberForm in Air (60s) b) FiberForm in N 2 (120s) c) 4” IsoQ initial Shape 

Figure 3 Arc Jet stagnation tests of FiberForm in Air and in Nitrogen. 


conditions. The stagnation heating and pressure in air were 440 W/cm and 49 kPa and in nitrogen the heating and 
pressure were 465 W/cm 2 and 47 kPa. 

The tests indicate that FiberForm does not significantly recess in the absence of oxygen such as a nitrogen 
heating environment, relative to the recession of FiberForm in the 20% oxygen environment of air (see Fig. 4). The 
measured heating rate (shown as a black square in Fig. 4b) agrees well with the prediction. 



a) Recession Distribution b) Heating Rate Distribution c) Shear stress Distribution 


Figure 4 FiberForm Recession in stagnation test with Air and Nitrogen. 


The good uniformity of recession, across the face of the test article, is attributable to the relatively uniform 
heating distribution across the face of the test article as the name IsoQ implies. It is interesting to note that the 
recession does not seem to be affected by the high rate of surface shear seen at the shoulder of the stagnation puck. 

The FIAT thermal response code, along with the FiberForm material properties database, does a good job of 
predicting the recession for both the air and nitrogen test cases. The small amount of recession in nitrogen is partly 


3 

American Institute of Aeronautics and Astronautics 



due to sublimation of the carbon and partly due to nitridation (N + C -> CN). These tests demonstrate that the 
primary recession mechanism in stagnation is due to oxidation. Similar results were obtained for PICA (but not 
shown). 

B. Shear Tests on a Swept Cylinder 

Tests were performed on FiberForm in the shear 
environment of a swept cylinder. These tests were 
performed in both Air and nitrogen in the AHF at the 
same arc heater conditions that the stagnation IsoQ 
models were tested. The swept cylinder consists of a 
304mm long water cooled copper cylindrical segment of 
radius 101mm with a pocket in the aft half of the fixture 
where material samples were mounted (see Fig. 5). 

The model’s axis was inclined 40° relative to the 
centerline of the arc jet so that the flow would impact the 
face of the cylinder at a 40° grazing angle. The 

upstream end of the cylinder was mounted 140 mm off- _ ^ , 

center from the jet centerline so that its nose was outside Fl S ure 5 Swe P* C y |,nder n,ounted “ the arc * et 
the flow. The heating conditions on the downstream end of the swept cylinder were measured by Gardon Gages that 
were mounted in a water-cooled copper plate that sat in the pocket of the swept cylinder. The resulting heating and 
pressure distributions on both a 30° and 40° inclined swept cylinder are shown for the runs with air and nitrogen (see 
Fig. 6). 





a) Heating Rate Distribution b) Surface Pressure Distribution 


Figure 6 Heating and pressure distribution along the pocket region of the swept cylinder. 


The heating and pressure are significantly higher for the 40° cases. The pressure is seen to diminish with 
distance along the cylinder - this is due to the divergence of the jet. The jet impacts the face of the cylinder at a 
more normal incidence angle at the upstream end of the cylinder than it does at the downstream end of the cylinder 
(where the flow impacts at more of a grazing angle). Figure 1 1 in the next section illustrates this. 

FiberForm was tested in the 40° angle configuration in the interest of testing at the highest possible heating, 
pressure, and shear (shear is similar for the two cases). The resulting recession rate is shown Fig 7. Similar to the 
stagnation tests, the FiberForm does not recede much in nitrogen compared to the recession rate seen in Air. Similar 
results were obtained for PICA (not shown). FIAT under-predicts the recession of FiberForm in Air by 50%. This 
under-prediction is similar to that seen for PICA (Fig. 1). 

The fact that the FiberForm recedes minimally in nitrogen is an indication that mechanical erosion does not play 
a significant role in the recession of non-oxidized FiberForm. However, in an oxidizing environment (air flow) the 
material near the surface may be less dense, and possibly more fragile, as the fibers thin due to oxidation. It is 
conceivable that mechanical erosion may still play a role in the recession of carbon char in its weakened state 
resulting from oxidation. These tests by themselves do not rule out the mechanical erosion theory. 
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C. Shear Tests on a Wedge 

Additional tests were performed in the 60 MW Interactive 
Heating Facility 7 (IHF) arc jet at NASA Ames Research Center in 
a wedge test fixture (see Fig. 8). Test conditions were the same 
used for the MSL tests of PICA (aka MSL 3 condition 8). This 
was the highest condition tested by MSL, with the arc heater 
running at 6000 amps and a total mass flow rate of 0.85 kg/s 
through the 6” (152mm) nozzle of the IHF. The corresponding 
stagnation heating and pressure on a 50mm radius Hemi-spherical 
probe was -2000 W/cm 2 and -150 kPa (respectively) at a 
location 76mm downstream of the exit of the nozzle. The heating 
on a 20° half-angle wedge ranged from 350 W/cm 2 at the 
upstream end of the wedge to 200 W/cm 2 at the down stream end 
of the wedge. 

The FiberForm recession in air is largest near the upstream 
end of the test article and significantly less at the aft end of the 
test article (see Fig. 8c). A FIAT calculation of recession (at 
x=75 mm) is under-predicting the recession. The measured 
recession is about 50% greater than the FIAT calculation. This 1.5 ratio of measured-to-predicted recession is the 
same as that seen in PICA tested at these same conditions. The ratio of measured-to-predicted recession being the 
same for both PICA and FiberForm is an indication that the missing physics (causing the mis-prediction) is likely to 
be the same for both. This is good news because we can use the simpler FiberForm test case to help test theories as 
to why FIAT is under-predicting the recession. In particular we will be able to more easily test the theory that the 
diffusion of oxygen is not proportional to the diffusion of heat in these tests. 



Figure 7 Recession of FiberForm tested in the 
Swept Cylinder in both Air and Nitrogen 




a) Wedge in IHF 


b) Post Test FiberForm 
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c) FiberForm Recession 


Figure 8 Arc jet tests of FiberForm in a wedge. 


III. Analysis 

Analysis using the DPLR 8 CFD code and the FIAT material response code were used to help understand the 
recession rates seen in these stagnation and shear tests. 

A. FIAT Material Response Code 

Prediction of recession is routinely performed with the FIAT 5 material response code in conjunction with a 
material properties database. The FIAT code assumes equilibrium chemistry and boundary layer theory to relate the 
recession rate to the surface heat transfer rate. 


1. Surface Energy Balance 

Material response codes (such as FIAT) solve the following surface energy balance equation (in one form or 
another) derived by Rindal, Kendall, et.al. 9,10 ’ 11 ’ 12 

P,«,c„ ( n r - *„) + + m s h s +p, U ,C„{X(z‘ - r„)F - B’h „ }+ - FosX - = 0 
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Z t * and Z iw * are the mass fractions of species i evaluated at the boundary layer edge and at the wall (respectively). 
This (4 th term) represent the transport of the chemical energy across the boundary layer associated with the species 
diffusion. Note that species gradients through the boundary layer are caused by reactions at the wall and thermal 
gradients. B’ is defined as 

m n + m o (nv) 

B'= g = — 

Pe U e Cm Pe U e 

The energy equation simplifies when we assume that C M = C H and further assume that species diffusion 

Z . h t w = h w and } Z . h t w = h e 

i i 

p e u e C„[H r -h w ]+rit c h c + m g h g ~(pv) w h w +a v q„ -FasJT* -q con< = 0 


FIAT solves the above surface energy equation in the course of a transient analysis. The heat transfer coefficient 
C H contains the effects of blowing, and is modeled as 

ln(l + 2 ZB') 


C =C 

^ Ho 


2AB' 


(i) 


where C Ho is the heat transfer coefficient for the case without blowing. C Ho is typically obtained from CFD with 
input from the tests. 


2. Steady State Energy Balance 

In the case of steady state ablation where the pyrolysis zone leading edge is propagating inward into the virgin 
material at the same rate as the ablation front is propagating inward into the char layer, then the conduction term is 

equal to fh c h c 4- th g h g — ( m c + m g )AH™ K allowing for further simplification (also assuming that the view 


factor F= 1 and ignore radiation coming to the surface from the outside), 


P e u e C H [H r - (1 + B')h w + B'AH™«]-aeX = 0 

One can think of the energy equation as the equation that determines the surface temperature of the ablator. 
Rearranging the steady state surface energy equation above, and using the definition of C H gives 

f h A tt29% k \ 


- Qhw 


v 


H-h 


It is interesting to note that for PICA and FiberForm both B ' and (h w -AH^ 98K ) /(H r - h w ) are typically much less 
than 1 and their product is of order 10 2 . Consequently, q hw ~<r£ w T w 4 is a fairly good estimate of the hot wall heating 


rate, or conversely, surface temperature T w ~( q hw /(<j£ w ) ) 


0.25 


3. Carbon Mass Loss for a Non-Pyrolyzing Ablator 

For cases where there is no pyrolysis gas, such as with 
FiberForm, Graphite or Reinforced Carbon-Carbon, the 
analysis further simplifies. 

A Multi-component Ablation Thermochemistry code 
(MAT) 13 was used to develop Bf tables for pure carbon 
ablators (with no pyrolysis gas) in an air environment. The 
results are plotted in Fig. 9. It is no accident that Bf is a 
constant of 0.175 over a large range of temperatures and 
pressures. The mass fraction of oxygen in air is 0.23 
which when multiplied by the ratio of molecular weights, 

MJM 0 , gives 0.23*M c JM o = 0.175. Equilibrium thermal 
chemistry (JANNAF database), from which MAT 
determines the species concentrations at the wall, says that 
virtually all oxygen (at the wall) will react with carbon to 
create CO when the gas at the wall is in a temperature 
range from 1000 to 3000K. In other words, given enough 
time (i.e., thermal equilibrium) all oxygen will combine Figure 9 B’c values for a pure Carbon Ablator 
with carbon atoms, irrespective of the heating rate in a as a function of Temperature and Pressure 
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range from 10-450W/cm 2 . Above 3000K carbon will also sublime in addition to oxidize, but that is not an issue for 
these tests since they were all run at a heating rate of 450W/cm 2 or less with temperatures between 220(XT W <2900K 
Other values of B c * were obtained for cases Argon mixed in with the air - the arc jet injects small amounts of the 
inert species Argon to protect the electrodes from oxidation. 

Rather than running a full FIAT simulation we will be able to model the recession rate as a constant times the 
heat transfer divided by enthalpy. To see this we can write the carbon mass loss rate as 

„ Df „ ln(l + 2AB' 

m c ~ Pe U e^H^ ~ Pe U e^Ho ~ * 


Here the standard blowing rate correction is used (Eq. 1.) with C Ho being the heat transfer coefficient associated 
with the no-blowing case and A- is a transpiration coefficient (blowing reduction parameter) assumed to be 0.5 for 
laminar flows. Also we use B c * = 0 . 1 75{1-C ar ) where C ^ is the mass fraction of argon. The heat transfer coefficient 
for the non-blowing case is calculated from either measured or CFD predicted wall heating and flow enthalpy, using 
p e UeCHo = qw/(Hr-h w ). Knowing the density of the carbon material, the recession rate can be determined 


S = 



ln(l + 2/15' ) ) 

P ablator {P r ~ ) 

< 2^ J 


ip 


ablator 


(2) 


FIAT requires p e u e C Ho and H r as inputs and performs a similar calculation inside the code. The above equation is 
used to evaluate the FIAT equivalent recession rate without actually resorting to running the FIAT code. This is the 
recession rate commonly reported in the next section that is referred to as a FIAT equivalent. 


B. Carbon Oxidation Model within DPLR CFD Code 

The DPLR CFD code was modified to include a surface carbon oxidation model as one of its boundary 
conditions. The modification allows oxygen from the flow field to react with carbon from the surface to generate 
CO gas ( O + C surface -> CO ). This is accomplished by adding sink and source terms at the wall in the species 
transport equations for O and CO (respectively). The terms are as follows 


RT . M . 

W o = P^J -TTT- and = ~^r W o 

V 2 nM „ M n 


The surface normal blowing velocity comes from the mass addition of carbon from the TPS and is given by 

v=cr ^lJKl 

° M V 2nM n 


In this model, the reaction rate has a user controlled reaction probability, y. A y of 0.9 is used for the calculations 
reported in this paper - however tests with other values of y indicate that the overall solutions are relatively 
insensitive to this parameter. This insensitivity is due to the fact that the reactions are limited by the rate at which 
oxygen arrives at the surface rather than kinetics (which is relatively fast). 

The combustion products (CO molecules) diffuse into the boundary layer creating a layer that gets in the way of 
oxygen diffusion to the wall. The balance between oxygen diffusion to the wall and CO diffusion away from the 
wall, simulated in the gas portion of the CFD code, dictates the rate of recession. This process is referred to as 
diffusion limited reaction rates. 

In the course of testing the carbon oxidation model, additional surface interaction processes were investigated, 
such as oxygen and nitrogen catalysis (0+0^0 2 and N+N->N 2 ) in conjunction with 0+C(s)->C0. Also the 
oxidation by 0 2 +2C(s)^2C0 was tested. The addition of the 0+0->0 2 reduced the recession somewhat, while the 
0 2 +2 C(s)-> 2C0 reaction increased the recession slightly. The two reactions together resulted in a recession rate 
that was similar to that of using only the 0+C(s)^C0 reaction. For the purpose of this work we incorporated all 4 
reactions. 


C. DPLR Solution Process 

All arc jet simulations are performed in two steps. The first step is to compute the nozzle flow. The second step 
is to solve the flow around the test article and its grid after interpolating the nozzle flow solution onto the outer 
boundary of the test article grid. 
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The nozzle flow solution 
encompasses the upstream 
plenum and the downstream test 
box, following the approach 
described by Gokcen et al. 14 
(see Fig. 10). The facility- 
measured mass flow rate and 
bulk enthalpy are imposed at 
the upstream end of the plenum. 

While the downstream test box 
pressure is prescribed at the 
walls of the test box based on 
the measured test box pressure 
during the time the model is in 

the flow. The degree to which Figure 11 Nozzle solution including upstream plenum and downstream test box. 

the jet blooms is a function of 

the box pressure. Fig. 10 shows the result of a simulation of the 
AHF with a swept cylinder model located in the jet. 

The test article simulation (swept cylinder in this example) is 
performed as a local problem, in which the nozzle solution is 
interpolated onto the outer edge of the test article’s flow domain. 

A close up of the cylinder within the nozzle flow domain is 
shown in Fig. 11. By interpolating the nozzle flow solution onto 
the test article domain we are replicating the varying radial 
distribution of the nozzle flow. The swept cylinder is particularly 
sensitive to the nozzle flow, since the cylinder’s nose sits outside 
the arc jet. This causes the cylinder to develop a stagnation point 
on the face of the cylinder with most of the flow traveling in the 
stream-wise direction, and a small portion of the flow traveling 
toward the nose of the cylinder. Coupling the nozzle flow 
solution to the swept cylinder solution is necessary to accurately 
simulate the flow over the swept cylinder model. 



Figure lOClose up of the swept cylinder 
domain within the nozzle domain. 



IV. Carbon Oxidation Model Simulations 

DPLR simulations, using the Carbon-Oxidation boundary condition are performed for the tests described in 
Section II. 

A. Stagnation Simulations 

In stagnation flows we expect to see good agreement between carbon-oxidation modeling and heat-transfer 
based modeling since the thermal boundary layer and the CO boundary layer start at the same location (the 
stagnation point) and grows at similar rates. Thus we can expect that the heat transfer based model (Eq. 2) to match 
the carbon oxidation model results. CFD simulations were performed for the stagnation tests of FiberForm in air 
(see Fig. 12). 

As expected, the FIAT equivalent recession rate is similar to that predicted by DPLR using the carbon-oxidation 
boundary condition. Both predictions of recession agree reasonably well with the measured recession rate. The 
FIAT equivalent calculation agrees nicely with the full FIAT prediction (shown as a blue square in Fig. 12d). The 
full FIAT calculation includes a 5% fail component to recession that is not included in the FIAT Equivalent 
calculation and thus the slight difference. 

CFD calculations of a CEV capsule simulation (not shown) also show that the FIAT equivalent and carbon 
oxidation model give the same prediction of recession. This is good news for heat shield designers who rely on 
FIAT and other similar codes to predict recession on stagnation geometries such as capsules. 
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B. Wedge Simulations 

The CFD simulation contains two zones; one shown in green (lower left comer of Fig. 13) which contains a cold 
isothermal copper boundary condition, and a second zone (in red) which has the FiberForm carbon-oxidation 
boundary condition. The FiberForm recession rate prediction using the DPLR carbon-oxidation model is shown in 
red in the lower right hand comer of Fig. 13. Also shown in the lower right hand comer (in green) is the equivalent 
of a FIAT calculation of recession (a heat transfer coefficient based prediction). The difference between the two 
predictions of recession is due to the difference between heat transfer coefficient and oxygen mass transport 
coefficient inherent in flows of this nature. The carbon-oxidation model is in much better agreement with the 
measured recession (shown as squares in Fig. 13) than is the FIAT equivalent prediction. 

Calculations were also performed on an all carbon wedge (not shown), in which the prediction of recession by 
the carbon-oxidation model compared very well with the FIAT equivalent prediction. This is further indication that 
stagnation flows (nose of wedge) are not an issue so long as there is no abrupt change in materials. 
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Figure 13 DPLR Prediction of FiberForm in a Wedge Configuration (no shape change). 


The remaining difference between the measured recession and the recession predicted by the carbon-oxidation 
model may be due to shape change effects. The thought is that the dished out shape of the recessed test article can 
create a low heating rate (and low oxidation rate) in the upstream region on the test article, and a high heating rate 
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Figure 15 DPLR Prediction of FiberForm in a Wedge 
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(and high oxidation rate) in the downstream 
region. This is due to expansion and 
compression waves generated by the test 
article’s cavity shape. To test this the DPLR 
code was further modified to allow for shape 
change as a function of time. At periodic 
intervals (every 1000 iterations after 
convergence to a steady state solution) DPLR 
was programmed to translate the FiberForm 
surface according to the predicted rate of 
oxidation multiplied by a time step, At. The 
grid was adapted (by stretching) to conform 
to the new geometry and the solution was re- 
converged. This process was repeated 21 
times with Ati=l second time intervals to 
achieve the final shape (after 21 seconds) 
shown in Fig. 14. In the implementation, a 
recession limiter was applied at the upstream 
end of the FiberForm test article to prevent 
the formation of a backward facing step 
which was incompatible with the grid 
topology - only applied to the first 5 grid 
surface grid points of the test article. Also 
there was no attempt to adapt the grid to the 
compression waves emanating from the 
downstream forward facing ramp. 

The final shape produces a lower pressure and heat transfer in the upstream region of the recessed shape where 
flow expansion takes place downstream of a backward facing ramp, and a slightly higher pressure and heat transfer 
in the downstream region of the recessed shape where compression takes place on the slightly inclined forward 
facing ramp. The recessed shape allows the oxygen rich air to skip over the cavity and impinge on the downstream 
region of the flow causing a relatively high recession rate downstream. The final recessed shape seen in Fig. 14d 
agrees well with the measurements in the 
downstream region of the flow and under- 
predicts the recession in the upstream 
region of the flow. 

It would appear that the non-flight like 
nature of the test article (cold copper 
upstream of carbon) is at least partially 
responsible for the high recession rate seen 
in tests relative to FIAT predictions. Most 
of the anomalous recession behavior seen 
in shear tests may simply be an artifact of 
the test configuration (copper upstream of 
test article). 

C. Swept Cylinder Simulations 

Similar to the wedge, the swept 
cylinder has a water-cooled copper section 
upstream of the test article (FiberForm in 
this case), so we might also expect to see 
some effect of the sudden jump in 
boundary condition (from cold copper to 
hot carbon). In this case (Fig. 15) the 
difference between the carbon-oxidation 
model and “FIAT equivalent” is less 
pronounced than it was with the wedge. 
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Figure 14 DPLR Prediction of FiberForm in a Swept Cylinder 
Configuration (without shape change) 
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The mismatch in recession rate between the FIAT equivalent and the carbon-oxidation model is greatest at the 
upstream end of the swept cylinder (as it was in the case of the wedge). After a short distance downstream of the 
copper-FiberForm junction the two predictions converge to nearly the same answer. This short distance to achieve 
similarity may be a result of the nature of the swept cylinder. That is to say that the upstream portion of the flow is 
more readily washed over the side and the centerline flow is replenished with flow from the free-stream. Another 
difference between the wedge and swept cylinder is the inclination angle (20° in the case of the wedge and 40° in the 
case of the swept cylinder. Because the swept cylinder is more normal to the flow than the wedge, the swept 
cylinder more closely resembles a stagnation flow with a large component of the flow avoiding contact with the 
copper section upstream. Never the less, the carbon oxidation model agrees slightly better with the data than does 
the “FIAT equivalent” prediction. 

Shape change effects on the swept cylinder were considered in Fig. 16. CFD predictions of heat transfer, 
pressure distribution and total recession are shown at time t=21 sec into a simulation with time varying shape due to 
recession. The recession rate predicted by the carbon-oxidation model at each time step is used to determine the 
increment of recession (during a 1 
second interval) that is applied to the 
surface location. 

After 21 seconds, the shape of the 
surface is sufficiently recessed to 
produce a significantly concave region 
behind the copper-carbon junction. 

The flow appears to separate at the lip 
of the backward facing ramp on the 
upstream end of the FiberForm. The 
heating and pressure is relatively low 
in this concave region. Downstream 
of the separation the flow reattaches, 
producing higher heating and pressure 
than seen in the case without recession 
(Fig. 16). Incorporating shape change 
into the carbon oxidation model 
simulation improves the agreement 
with the data somewhat, but not 
enough to explain all of the difference 
between the prediction and the data. 

V. Conclusions 

. This paper describes an effort to understand the relatively high rate of recession seen in shear tests relative to 
FIAT predictions. To this end, arc jet tests were performed on FiberForm, which is a light-weight carbon preform 
used as a precursor to the fabrication of PICA. The tests showed that FiberForm recession was similar to that of 
PICA when adjusted for the difference in density. Similar to PICA, the measured recession of FiberForm is higher 
in a shear environment than FIAT would predict. Several theories, attempting to explain the high recession rate, are 
explored in this paper. The data is proving useful for understanding the high recession rates seen on PICA in shear 
tests. 

One theory, mechanical erosion, was investigated in a pure nitrogen arc jet test. The measured recession rate of 
FiberForm in pure nitrogen arc jet tests (in both stagnation and shear) was negligible compared to similar tests in air. 
This is an indication that the non-oxidized char layer of FiberForm (nearly identical to surface char layer of PICA 6 ) 
is not fragile to the point that the shear forces can remove material by mechanical erosion (at the levels tested here). 
It is conceivable that mechanical erosion may still play a role in the recession of carbon char that is weakened by 
oxidation, but we think this is unlikely. While we can not rule out the mechanical erosion theory, these tests are a 
necessary step toward doing so. 

A second theory, that the diffusion of oxygen is not proportional to diffusion of heat was tested analytically. The 
theory being that the non-flight like test configuration involving a copper test fixture upstream of the carbon test 
article is responsible for the high rate of recession. In particular, the copper fixture produces an enthalpy depleted 
boundary layer, while leaving the boundary layer oxygen rich. The dissimilar profiles of enthalpy and oxygen 
concentration produce a condition that violates the assumption in FIAT that C m = C h . Shear test simulations with 
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Figure 16 DPLR Prediction of FiberForm in a Swept Cylinder 
Configuration (with shape change) 
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DPLR and a carbon-oxidation model produced better predictions of FiberForm recession than did the FIAT 
equivalent prediction. It would appear that the high rate of recession seen in shear flows can at least partially be 
explained by the non-flight like test configuration involving a copper test fixture upstream of the carbon test article. 
As a cross-check of the carbon-oxidation model, stagnation simulations were performed, showing that in stagnation 
flows there is no significant difference between carbon-oxidation model and the FIAT equivalent. 

A third theory, test-article shape change, was investigated. The thought being that the dished out shape of the 
recessed test article can create a low heating rate (low oxidation rate) in the upstream region on the test article, and a 
high heating rate (high oxidation rate) downstream region. This is due to expansion and compression waves 
generated by the test article’s cavity shape. CFD simulations, which included time varying shape change of the test 
article were performed and confirmed the aforementioned theory, but the changes in heating (and predicted 
recession rate) were not large. The greatest effect of shape change was in the upstream region where the flow was 
inclined to separate from the upstream lip of the dished out shape. As a result of the oxygen rich flow was able to 
skip over this region and impinge on a region further downstream, slightly increasing the recession rate downstream 
relative to the recession rate of the original flat shape. 

Effects that may be important, that were not modeled in this study include material porosity and Gortler vortices. 
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